Time-reversal symmetry breaking in circuit- QED based photon lattices 
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Breaking time-reversal symmetry is a prerequisite for accessing certain interesting many-body states such as 
fractional quantum Hall states. For polaritons, charge neutrality prevents magnetic fields from providing a di- 
rect symmetry breaking mechanism and similar to the situation in ultracold atomic gases, an effective magnetic 
field has to be synthesized. We show that in the circuit QED architecture, this can be achieved by inserting 
simple superconducting circuits into the resonator junctions. In the presence of such coupling elements, con- 
stant parallel magnetic and electric fields suffice to break time-reversal symmetry. We support these theoretical 
predictions with numerical simulations for realistic sample parameters, specify general conditions under which 
time-reversal is broken, and discuss the application to chiral Fock state transfer, an on-chip circulator, and 
tunable band structure for the Kagome lattice. 

PACS numbers: 42.50.Dv, 42.50.Ct, 71.36.+C 
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I. INTRODUCTION 

Since the first pioneering papers in 2006 d-d], theoretical 
interest in the many body-physics of interacting photons or 
polaritons in lattices has flourished. Such photon lattices, see 
Fig.Q]for an example, are perceived as an interesting venue for 
quantum simulation J3l and for studying strongly correlated 
systems composed of polaritons yrtZfl. Hopes are that, once 
realized in experiments, such systems could complement the 
achievements in research with ultracold atomic gases 
which are currently leading the charge. 

Much recent work has focused on the quantum phase transi- 
tion between polaritonic Mott-insulating and superfluid states 
using various approaches j1f-{3l [l0l - [r7tl . and at this point there 
seems little doubt that the quantum phase transition is in the 
same universality class as its counterpart in the Bose-Hubbard 
model lfl8l - l20Tl . It is thus natural to ask, what physics beyond 
Bose-Hubbard might photon lattices have to offer? 

Recent work by several groups has highlighted the interest- 
ing implications of dissipation and external driving, and thus 
promoted the quantum phase transition to a nonequilibrium 
phase transition between different possible steady states I2TI - 
12411 . A second route to physics beyond Bose-Hubbard, is to 
explore phases with broken time-reversal symmetry, of which 
fractional Quantum Hall phases are the most celebrated exam- 
ple jum. 

To access such phases, a technique for breaking time- 
reversal symmetry is required. In contrast to electron gases, 
but similar to ultracold atomic gases l27T - [30ll . polariton sys- 
tems face a challenge when trying to break time-reversal sym- 
metry: due to the charge neutrality of polaritons, an exter- 
nal magnetic field cannot readily be used to achieve break- 
ing of time-reversal, and instead an effective magnetic field 
has to be synthesized. A first proposal for cavity arrays with 
trapped three-level atoms and involving ac driving with spe- 
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FIG. 1: (Color online) The Jaynes-Cummings lattice as an exam- 
ple of a photon lattice. Its circuit QED realization would consist 
of superconducting resonators (e.g., coplanar waveguides, schemati- 
cally shown as rectangular boxes), each of which would be coupled 
to a superconducting qubit (symbolized as dots centered in the res- 
onators). Microwave photons would hop between nearest-neighbor 
resonators, with the coupling strength k set by the mutual capaci- 
tance between resonator ends. Interaction between the photons and 
the superconducting qubits with strength g would induce an effective 
photon-photon interaction. 
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cific phases has been published by Cho et al. 13 1H . In addition, 
photonic edge states and analogs of the quantum Hall effect in 
photonic cr ysta ls have recently been investigated by Haldane 
and Raghu ll32l l33ll and also probed experimentally lllill . 

In the present pap er, we demonstrate that in the circuit QED 
architecture 035143711 breaking of time-reversal symmetry can 
be achieved by inserting simple superconducting circuits into 
resonator junctions and applying purely dc electric and mag- 
netic fields. In our scheme, photons are transferred from res- 
onator to resonator via virtual intermediate excitations of cou- 
pler circuits. We expect that the use of passive coupling ele- 
ments and the absence of any ac fields pumping internal lev- 
els may avoid some of the challenges posed by dissipation. 
Our analysis shows that for broken particle-hole symmetry 
(caused by a dc electric field), polaritons can acquire an ef- 
fective gauge charge and hence become susceptible to an ex- 
ternal magnetic field so that time-reversal symmetry is bro- 



2 




FIG. 2: Basic scheme of a three-port coupling element, connected 
capacitively to three transmission-line resonators with annihilation 
operators ay for photons in the relevant mode of the resonators enu- 
merated by j = 1, 2, 3. 



ken. We emphasize that such passive coupling elements cor- 
respond to an important step towards substituting commercial 
microwave circulators with on-chip circulators much smaller 
in size. This could pave the way for integrating circulators 
into larger arrays of resonators and could open interesting and 
new perspectives for correlated polariton systems. 

The remainder of the paper is organized as follows. In Sec- 
tion|II]we explain the generic consequences of integrating pas- 
sive coupling elements into a resonator array and using them 
to break time-reversal symmetry. The passivity condition al- 
lows us to adiabatically eliminate the coupling elements and to 
obtain an effective photonic tight-binding model with broken 
time-reversal symmetry. We emphasize the gauge-invariant 
phase sum (mimicking the contour integral of the magnetic 
vector potential in the continuous case) as a useful concept 
for determining whether time-reversal invariance holds. Ap- 
plications of such coupling elements, including the prospect 
of an on-chip circulator conclude the section. 

Section HH1 then details our proposal for a physical realiza- 
tion of passive coupling elements in the circuit QED architec- 
ture. Specifically, we consider a system consisting of coplanar 
waveguide resonators which capacitively couple to small su- 
perconducting rings interrupted by three Josephson junctions 
("Josephson rings"), which are inserted into the junctions be- 
tween resonators. Using circuit quantization, we derive the 
Hamiltonian of this system and discuss the diagonalization of 
the Josephson rings. 

In Section HVl we finally show how the adiabatic elimina- 
tion of the ring degrees of freedom yields an effective photon 
Hamiltonian of the desired type. We discuss the general re- 
quirements for achieving time-reversal symmetry breaking in 
this scheme, and present results from numerical simulations 
which underline the proposal's feasibility with realistic device 
parameters. 

We end with conclusions and an outlook in Section [V] 
Some additional details of calculations and a self-contained 
summary of time-reversal symmetry in quantum mechanics 
are provided in several appendices. 



II. PASSIVE COUPLING ELEMENTS FOR BREAKING 
TIME-REVERSAL SYMMETRY 

For the general discussion of breaking time-reversal sym- 
metry by utilizing virtual excitations of a coupler circuit, we 
consider a junction composed of three resonators ffda coupled 
to a central "circulator" system, see Fig. [2 and described by a 
generic Hamiltonian of the form 



H B . 



(1) 



Here, a, and aj (j = 1,2,3) are annihilation and creation op- 
erators for photons in the relevant mode of resonator j, with 
corresponding (angular) frequency ui r , (Note that throughout 
the paper we use units with h = 1.) The capacitive coupling 
between resonators and the degrees of freedom Bj of the cou- 
pling element is described by the second term in Eq. (Q]). 

We shall assume that the coupling element remains pas- 
sive, i.e. the coupler only transfers photons via intermediate 
virtual excitations and otherwise remains in its ground state at 
all times. Consequently, the coupler degrees of freedom can 
be integrated out (or, in other words, eliminated by a canon- 
ical transformation of Schrieffer- Wolff type |H,[39j]) so that 
one obtains an effective photon Hamiltonian _ff e ff {dj , The 
details of the effective Hamiltonian H e g generally depend on 
the specific realization of the passive coupling element, and 
we will go through the explicit derivation of H e s for the cir- 
cuit QED realization we propose in Section[III] Here, we first 
explore the generic properties of the effective photon Hamil- 
tonian. 

We are interested in a passive coupling element that does 
not destroy the three-fold symmetry of the system. As a re- 
sult, there is a gauge in which H e s is invariant with respect 
to cyclic permutations of the indices j = 1, 2, 3. Further, we 
assume that H e g allows for hopping of photons between res- 
onators, but does not induce photon-photon interaction. (This 
assumption is realistic, as we show in SectionHIIl) As a result, 
H e ff is anticipated to be a quadratic form of the annihilation 
and creation operators a,j, a]. Explicitly, the Hamiltonian will 
take the form 



H. 



eff 



a 3 a 2 



a 2 a\) + H.c. 



E 

3=1 



oj r a' J aj, (2) 



where uj' t denotes the resonator frequency (possibly includ- 
ing a renormalization), and t = Ke %v (k = \t\ > 0) is the 
complex- valued hopping matrix element for photons 17711 . 

When does the effective Hamiltonian (|2]) describe the situa- 
tion of broken time-reversal symmetry and when does time- 
reversal symmetry remain intact? Formally, time-reversal 
symmetry holds whenever the time-reversal operator leaves 
the Hamiltonian invariant, i.e. QHQ^ 1 = H lf78ll . As de- 
tailed in AppendixfA] for the present case this is true if there 
is a gauge transformation of the form 



(3) 
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which makes the Hamiltonian real-valued when represented 
in the photon number basis. For the three-resonator junction, 
the existence of such a gauge transformation is checked as 
follows. According to Eqs. (O and ([3J, an attempt to find 
a gauge transformation to make the Hamiltonian real-valued 
leads to the three equations 
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where z±, Z2, 23 € Z are arbitrary integers. These equations 
for the gauge phases ipi, ip 2 and can only be solved (and 
hence time-reversal symmetry is intact) if the condition 

3cp = zn (ze Z), (5) 

obtained by summing the three equations (01, holds. Thus, for 
the present case of a three-resonator junction we find: time- 
reversal symmetry is intact if and only if ip <E ^Z. 

To extend this statement to general photon lattices with 
more resonators, 

iZe ffi = ^2 '<.;"<", + X! Urd^cij (tji = t*j) (6) 

3 

it is important to identify the phase in Eq. © as a gauge- 
invariant quantity, which for discrete lattices plays a role anal- 
ogous to the contour integral § ds ■ A of the vector potential A 
in the continuous case. [For simpler notation the prime in uj' r 
has been dropped.] We write the gauge-invariant phase sum 
in the form 

^2 <Pij ■= arg Y[ tij, (7) 
CM cm 

where C specifies a closed path in the discrete lattice; see Fig. 
[3] for an illustration. In these terms, the statement of Eq. (© 
can be extended to larger systems where time-reversal sym- 
metry can be shown to be intact if and only if the gauge- 
invariant phase sum is an integer multiple of ir, 

ipij e 7rZ (8) 

CM 

for any closed lattice path C. 

To illustrate the implications of broken time-reversal sym- 
metry, we discuss three examples: the clockwise or coun- 
terclockwise state transfer of a single photon Fock state be- 
tween resonators, circulator behavior for signals propagat- 
ing in semi-infinite transmission lines, and tunability of the 
Kagome tight-binding band structure. These examples are re- 
alizations of the simplest setting possible: resonators coupled 
via coupling elements and without any photon-photon interac- 
tion. The fascinating scenario of systems of interacting pho- 
tons with broken time-reversal symmetry is beyond the scope 
of this article and will be addressed in a future paper. 



FIG. 3: (Color online) Illustration of the gauge-invariant phase sum 
around a loop, W,c{ij] W = V12 + ¥>23 + ^34 + <p4,u here for a 
particular plaquette C in a two-dimensional quadratic lattice. 

A. Chiral transfer of photon Fock states 

We consider the 3-resonator junction depicted in Fig. [2] and 
described by the effective Hamiltonian H e s, Eq. (f2]i. can 
be understood as a miniature tight-binding model with peri- 
odic boundary conditions. The eigenstates of H e ff are gener- 
ated by the creation operators 

4 = ^E^ /3 a] (9) 

and have corresponding eigenenergies 

Q k = w r + 2KCOs(27rfc/3 + ^). (10) 

Here, 27rfc/3 (k = —1,0, 1) are the allowed wave numbers 
in the first Brillouin zone. Recalling from Eq. © that time- 
reversal symmetry only holds as long as p e |-Z, it is not 
surprising that the simplest case of broken time-reversal sym- 
metry (where the energy spectrum set by 0^ becomes equidis- 
tant) is realized when p = ±7r/6, i.e. halfway in between the 
time-reversal symmetric points ip = and ±7r/3. 

To understand the effect of broken time-reversal symme- 
try, let us consider the dynamics of the system inside the one- 
photon subspace. We initialize the system in a Fock state with 
a single photon inside one resonator, say resonator j = 1, 
and follow its subsequent evolution in time. The evolution is 
obtained by solving the time-dependent Schrodinger equation 
with initial condition | ip(t = 0) ) = a] =1 | ). By using the 
inverse of the discrete Fourier transform in Eq. (0, the evolu- 
tion for ip = ±7r/6 is readily found to be 

1 1 

j ^(t) ) = — e iUrt e jfcv%tt - 2 ™ fc / 3 I Vfc ) , (ID 
v 3 fe =_i 

where \tpk) = A\. | ) denotes the single-photon eigenstates 
of iJ e ff- The dynamics may be visualized by plotting the prob- 
abilities 

p,(t)H(0K =1 |v<W)l 2 (12) 

for finding the photon in resonator j, see Fig. [4] As ex- 
pected from Eq. ( fTTT ). the dynamics is periodic with period 
t = 27t/v / 3k. More importantly however, the breaking 
of time-reversal symmetry results in chirality: the photon is 
transferred from resonator to resonator either clockwise or 
counter-clockwise depending on the sign of ip = ±7r/6. 



4 



Y Y T 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



: \ / \ / 


i3 /* 












A 






-i — i ^-r— i 


/ \ 



t 

x1 



FIG. 4: (Color online) Time evolution of a single-photon Fock state 
in the presence of a coupler with phase ip — n/6. The quantum state 
at the initial time t — is a Fock state with one photon in resonator 
j = 1, and both resonators 2 and 3 in the vacuum state. The photon 
occupation probabilities Pj are plotted as a function of time and show 
how the photon is transferred around the loop in a direction specified 
by the sign in tp — ±7r/6. The evolution is periodic with period 
r = 27r/v3« and the initial state is transferred into Fock state of 
resonators 2 and 3 at times t = r/3 and t = 2r/3, respectively. 



B. On-chip Circulator 



Circulators are lossless microwave elements with three (or 
more) ports, and have the crucial property that a signal en- 
tering port j is fully transferred clockwise toport j + 1 (or, 
alternatively counter-clockwise to port 114011 . This behav- 
ior must involve breaking of time-reversal symmetry, which is 
typically accomplished by embedding magnetic material, e.g. 
ferrite, in the device. Commercial ferrite circulators are typ- 
ically large (> 1 cm) and their size would make it rather dif- 
ficult to include large numbers in a photon lattice. It is thus 
interesting to explore the design of an on-chip circulator, suf- 
ficiently small in size and easy to fabricate, such that it could 
be included in large numbers. In addition to being essential 
for breaking time-reversal symmetry in polariton lattices, such 
devices would find great practical application in the circuit 
QED architecture for quantum information processing. 



Hamiltonian 

3 



Ke' v (aial + + a^a\) + H.c. 



H =ui r a^aj 

3 

+ E E u i b )q b 3i - * E E \fi h 3i a ) - H -' 



3=1 9 



3=1 9 



(13) 



where bj q are the annihilation operators for the three transmis- 
sion lines j = 1,2,3, and q is the mode index. 

We divide the full Hamiltonian H = H e s + H t \ + H{ n 
into the effective photon Hamiltonian previously discussed, 
the contribution from the semi-infinite transmission lines, and 
the interaction between them. Next, we employ the diagonal- 
ization of 77 e ff, see Eqs. (O and ( [Tol l, and rewrite the coupling 
Hamiltonian H mt in terms of the eigenmodes Ak, 



1 3 1 

EE E 



V3 



g j=lfc=-l 



f qe -^ijk/3 bjqA f k _ H c 



(14) 

To calculate ingoing and outgoing fields, we use input-output 
theory ll4lll42ll . As usual, formal integration of the Heisenberg 
equation of motion for bj q , 



lUJ a 



1 1 

b jq + — f* £ e 2 ^A kl (15) 



k=-l 



yields solutions which can refer to either an initial state at time 
ti = to in the distant past, or to a final state at time ij = t\ in 
the distant future: 

b jq (t) =e- iu ^- u) b jq (ti) 

i r* 1 

+ / dr e - l ^ {t - T) f* V e 2 ^ k/3 A k {r). 
V3 J ti fc= _ x 

(16) 

Proceeding with standard input-output theory, we approxi- 
mate the coupling matrix elements f q as constants within the 
relevant frequency range near f2/-, and employ the Markov ap- 
proximation l4lll . We then plug Eq. (fT6b into the equation of 
motion for A k , 



A k = iQ k A k - ^ E E h^ %3k/ % q 



(17) 



q j=l 



Let us demonstrate that circulator behavior in the sense 
of microwave engineering can indeed be achieved with the 
model Hamiltonian H e s, Eq. (|2]). The actual physical realiza- 
tion within the circuit QED architecture will be discussed in 
Section [ill] For simplicity, we consider a setting where mi- 
crowave radiation is fed into the system by capacitively cou- 
pling semi-infinite transmission lines to the three resonators 
shown in Fig. [2] The full system is then captured by the 



and identify the input and output modes as 

V^P q 

Here, p is the transmission line density of states, and k' = 
2ir j /| 2 p defines the effective photon decay rate. Applying 
the Markov approximation to the remaining time integral 14211 . 
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FIG. 5: Circulator behavior. The plot shows the normalized outgo- 
ing power |&° ut /&i n | 2 for the three ports j = 1, 2, 3 under coherent 
driving of port 1 with frequency Ud when the phase <p of the coupler 
element is adjusted to 7r/6. For drive frequencies close to the res- 
onator frequency, the signal is transferred from port 1 to port 2. The 
bandwidth of this circulator behavior is set by the photon hopping 
rates (n = «'/2 = O.lav). 

one obtains 

3 

A k = - itl k A k - v^V3E e-^ ijk / 3 bf 

3=1 

fe'=-i j=i 

Analogous expressions for A k (but with a crucial sign change 
in the last term) can be obtained when substituting the outgo- 
ing fields. By subtracting from Eq. ( fTSI l the equations obtained 
when using the outgoing field in either port 1, 2, or 3, one can 
derive the following relation between ingoing and outgoing 
modes: 

l 

b T = 5 j n + V«73 e 27ryfe/3 A fc . (19) 
fe=-i 

Finally, we eliminate the dependence on the circulator modes 
by substituting the solutions to Eq. ( fTSI ), which in frequency 
space can be expressed as 

MM = , ffi* ,._ £ e-^bfH (20) 
i(u) - il k ) - k'/2 f-f J 

3 = 1 

In total, one thus obtains the relation 

bf[u]=^[u] (21) 

K l 1 3 e 2xi(j-j')k/3 

ij — i 

between the ingoing and outgoing fields. For coherent driving 
with frequency w^, the ingoing and outgoing fields are charac- 
terized by c-numbers (6™' out [a;,J) and the normalized outgoing 




FIG. 6: (a) Using three-resonator junctions one obtains a photon lat- 
tice with uniform hopping, and the resonators (depicted as rectan- 
gles) form a regular honeycomb pattern, (b) The corresponding pho- 
ton lattice is the Kagome lattice, a hexagonal Bravais lattice (primi- 
tive vectors Ai, A2) with three atoms A, B, and C in the primitive 
unit cell (parallelogram shaded in gray). Adding coupler circuits in 
the junctions breaks time-reversal symmetry and introduces a phase 
factor e ±lv in the photon hopping elements, where the sign depends 
on whether photons are transferred with or against the sense of rota- 
tion (circular arrows). 

power can be calculated from Eq. (f2TT >. Assuming a drive at 

only one of the input ports, say port 1, the normalized outgo- 

1 1 2 

ing power on port j is | (b° u[ {uid\) / (b'lluid]) \ ■ Note that results 
for driving on any other port can be obtained by cyclic permu- 
tation of the port indices) 

As shown in Fig. [5] the device shows clear circulator be- 
havior when choosing <p = ir/6 for the photon hopping phase. 
The circulator behavior is strongest when the drive frequency 
ojd is close to the frequency of the resonators ui r . The band- 
width of circulator behavior is set by k in the configuration 
considered here. The condition k = k'/2 is required to 
achieve 100% transmission, and zero reflection at the input 
port. 



C. Tlinable band structure 

Incorporating coupler circuits into larger arrays of res- 
onators is useful for several reasons. As mentioned before, it 
may provide access to strongly correlated states of interacting 
photons with broken time-reversal symmetry. However, the 
usefulness of coupler circuits, is not limited to the interacting 
case. When leaving time-reversal symmetry intact, coupler 
circuits enables one to vary the (real-valued) photon hopping 
strength in situ and thus to systematically explore the quan- 
tum phase transition between a photonic superfluid and Mott 
insulator ifH-kil [Iol - [l7ll . Finally, when breaking time-reversal 
symmetry both magnitudes and phases of the photon hopping 
elements become tunable, which can make the photonic band 
structure tunable as we will show now. 

We consider a two-dimensional resonator array with uni- 
form photon hopping strength. With the circuit QED real- 
ization in mind (see Section [Til] for details), resonators may 
be imagined as coplanar waveguides and uniform coupling is 
readily achieved by using junctions composed of three res- 
onators at 120° angles. In this case, the coplanar waveguide 
resonators form a regular honeycomb pattern as shown in Fig. 
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13 a). Each resonator, depicted as a rectangle, represents a sin- 
gle lattice site. Thus, marking the center of each resonator as 
a lattice site and connecting nearest neighbor sites, one finds 
that the photon lattice is a Kagome lattice ll43ll . see Fig.[6jb). 

We briefly note that, due to its novel properties and physical 
realizations, the Kagome lattice has played an important role 
in various contexts of strongly correlated systems and frus- 
trated spin systems. Ferromagnetic and anti-ferromagnetic 
Ising 144144611 and Heisenberg II47I449I1 models have been stud- 
ied on the Kagome lattice. For the Hubbard model, the 
Kagome lattice is known to lead to flat-band magnetism ll50l - 
l52tl . The possibility to create optical Kagome lattices has 
also created interest in exploring this physics with ultracold 
atoms lf53ll54ll . Even more recently, the (fermionic) Hubbard 
model on the Kagome lattice has been revisited and shown to 
give rise to interaction-induced topological phases ll55l l56ll . 



Here, we show that even in the absence of interactions, the 
Kagome lattice displays an interesting tunable band structure 
when time-reversal symmetry is broken. The lattice particles 
(for us, photons) can assume eigenstates localized on only a 
few sites, giving rise to flat bands. Tuning the phase of the 
photon hopping makes it possible to modify the Kagome band 
structure and to switch the flat band to the top, middle or bot- 
tom band at will. 

To demonstrate this, we consider the tight-binding model 
of the Kagome lattice with nearest-neighbor coupling. The 
Kagome lattice is generated by a hexagonal Bravais lattice 
with primitive vectors Ai = a(l,0) and A 2 = %(l,y/S). 
The primitive cell contains three sites located at ro =0 (A), 
ri = Ai/2 (B), and r 2 = A2/2 (C), where positions are 
expressed relative to the origin of the primitive cell. The cor- 
responding tight-binding Hamiltonian is 



(22) 



^ y ] Cn m A nm + B^mCnm + A\ im B nm + C nm _yA nm + B*l l _ 1 _ m+1 C nm + A^ n+lm B n 



He, 



where we have already accounted for the fact that coupler circuits may introduce photon hopping with a complex phase factor, 
t = \t \ e tlf '. Working in reciprocal space, we find that the dispersion e s (k) of the three bands s = 1,2,3 is obtained from 
diagonalization of the following 3 x 3-matrix: 

/ u 2i*cos(k- Ai/2) 2tcos(k- A 2 /2) \ / A k \ 

H = J2(At Bl Ct)\ 2tcos(k-A 1 /2) u 2t* cos[k ■ (A 2 - Ai)/2] B k . (23) 

k \ 2t* cos(k- A 2 /2) 2tcos[k- (Ai - A 2 )/2] oj J \ C k J 

Compact analytical expressions for the band structure can be obtained for the special values tp = 0, tp = ir/6 and tp = tt/3: 

tp = 0: ei(k)=w-2*, e 2 , 3 (k) =u + t±\t\ + 8cos[ik • Ai] cos[|k • A 2 ] cos[|k • (Ai - A 2 )], (24) 
<p = ir/6: e 2 (k)=w, e lj3 (k) =u±2\t\Jl + 2cos[ik • Ai] cos[ik ■ A 2 ] cos[ik • (Ai - A 2 )]. (25) 



The case <p = tt/3 can be show to be equivalent to tp = it and 
is obtained from Eq. (f24t by switching the sign of t. The band 
structure for tp S is thus familiar from previous work, see 
e.g. Ref. lf57ll . Some of the results on the tight-banding band 
structure with broken time-reversal symmetry and an evalu- 
ation of the bands' Chern numbers have also recently been 
published in JH- 

The Kagome band structure for zero and nonzero <p is de- 
picted in Fig. 13 [For a beautiful discussion of the Dirac 
points in the band structure, visible in Fig. [7fc), and the lifting 
of the degeneracy by breaking time-reversal symmetry, see 
Ref. i32tl .l The characteristic flat bands occur exactly when 
ip E and depending on the specific phase, the flat band 
takes the role of the bottom or top band {tp = and tp = 7r/3), 
or that of the middle band (ip = n/Q and tp = ir/2). We note 
that phase values tp ^ [0, 27r/3) can always be mapped back 



into this interval via gauge transformations. 

Band flatness and the corresponding zero group velocit y ar e 
directly related to the existence of localized states ll52l l57ll . 
First, consider the phase values tp = ir/Q,ir/2 where the mid- 
dle band is flat. For periodic boundary conditions with a total 
of N primitive cells, the flat band corresponds to N energy- 
degenerate states. This degenerate subspace is spanned by the 
localized hexagon states \ip n ), where \ip n ) i s defined as the 
eigenstate localized on the n-th hexagon in the Kagome lattice 
with wavefunction amplitudes 

(jn|VO =(-l)V^/ 3 (26) 

on the six consecutive sites j = 0, 1, . . . , 5 of the hexagon. 
(Note that the | tp n ) states are linearly independent but non- 
orthogonal.) 

When the flat band is the top (or bottom) band, the situation 
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FIG. 7: Band structure of the Kagome lattice with complex hopping elements t — \t\ e' v for (a) ip = 7r/6, (b) <p = 7r/4, and (c) p = ir/3. 
In the top panels, the dispersion (e s — uj) of the three bands s = 1, 2, 3 is plotted in units of \t\. The first Brillouin zone corresponds to the 
hexagon centered at k = 0. The bottom panels show cuts of the dispersion along axes of high symmetry (see inset). For phases <p £ -|Z, the 
band structure exhibits flat bands. The position of the flat band can be switched from (a) middle to (c) top to bottom [for p = 0, obtained from 
(c) by reflecting all bands at (e — uj) — 0] by varying the phase <p. 



is slightly more complicated since the flat band touches the 
middle band at the k = point and the degenerate subspace 
is (N + l)-dimensional. The localized hexagon states with 
amplitudes (jn \ ip n ) = (— 1) J are eigenstates, but are not lin- 
early independent since their k = superposition J2 n I A» ) 
is identically zero. The localized states can be shown to span 
an (N — 1) -dimensional subspace, and the missing two k = 
states are obtained as 

| k = 0; 1 ) = -±= (4« - B ln) I ) , (27) 

* ran 

| k = 0; 2 ) = -±= Y,(Al n O | ) . (28) 

* ran 

The existence of localized photon states and the tunability of 
its band structure make the Kagome lattice with variable phase 
factors an interesting system for future experiments. Fur- 
ther theoretical studies will address the interesting question 
of strongly correlated states induced by photon interactions, 
which are expected to be non-perturbative in the presence of 
the flat band degeneracies of the Kagome lattice. 

III. PHYSICAL REALIZATION IN THE CIRCUIT-QED 
ARCHITECTURE 

Following the general discussion of broken time-reversal 
symmetry in photon lattices, we now turn to a concrete pro- 
posal on how to realize this physics in the circuit-QED archi- 
tecture. The essential idea is to insert superconducting cir- 
cuits into the junctions between resonators. These circuits 
then serve as coupling elements that transfer photons from one 
resonator to another and may break time-reversal symmetry. 

Our analysis will be organized into three subsections. In the 
first one, Section IIII Al we present the appropriate tools for 



index meaning 

i 6 {1, • . • , N} index decomposing resonator into LC elements 

j £ Z Josephson ring index 

k £ N excitation index for Josephson ring 

A £ Z resonator index 

jti, fixj £ {1,2, 3} component of ring j coupling to resonator A 

v £ N resonator mode index 



TABLE I: Summary of conventions for indices and their meanings, 



as used throughout Sections ITTIland lTVl 



modeling a transmission-line resonator capacitively coupled 
to arbitrary circuits at its two ends. We show how to system- 
atically obtain the exact eigenmodes of the resonator when it 
is coupled to arbitrary circuits at its two ends. These exact 
eigenmodes are then utilized in Subsection lHIBI to obtain the 
full Hamiltonian of a resonator array including coupling cir- 
cuits. Circuit quantization ll59ll allows one to switch to the 
quantum mechanical description of the full system. 

Notation in this section is heavy due to different types of 
objects (resonators, Josephson rings, etc.) that need to be enu- 
merated, and we have made every effort to be consistent in 
our naming of indices. For reference, the different labels are 
summarized in Table U 



A. Exact resonator eigenmodes in the presence of coupling 

We consider a system consisting of a transmission line cou- 
pled capacitively at its two ends to circuits described by La- 
grangians C' L R . The general configuration is depicted in Fig. 
[8] The Lagrangian of the full system can be cast into the form 
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N N 
AT 



(£g + \C a 4> 2 a ) -C L 4>i4>L - Cr^n^>r + \Cl4>x + \^r4> 2 n + — c<i 



i-lj 

1 w 

51 



(29) 



a=L.R 



i=l 



2ldz 



i=2 



Cl + Cr 



where the contributions Cl.r describe the circuits to the left 
and right (now including an additional capacitive contribution 
~ Cl,r due to the coupling), and C t \ the transmission-line 
resonator, modeled by an array of LC oscillators with capac- 
itances c dz and inductances I dz, where c, I denote the ca- 
pacitance and inductance per unit length. The capacitive in- 
teraction between resonator and attached circuits is denoted 

£int- 

It is useful to rewrite the transmission-line Lagrangian in 
compact matrix notation, 

U = ^ Ttf. - ] ^ 'y<f>. (30) 

with0 T = (0i, . . .,(/)n), 

(T)u> = 5u> {cdz + C L S n + C R 5 iN ) , (31) 



d e 




FIG. 8: (a) Transmission-line resonator attached through capacitors 
Cl and Cr at the left and right ends to arbitrary circuits. Panel 
(b) shows the dissection of the transmission line (capacitance and 
inductance per unit length denoted by c and t) into a series of LC 
circuits. The generalized flux variables adjacent to the resonator are 
given by <j)L and 4>r. 



and 



V 



1 

Tdz 



( 1 -1 

-1 2 -1 

-1 2 



V 



-1 2 -1 
-1 1/ 



(32) 



Generally, the eigenmodes <j> = <^,a„e _i " ,yt of the trans- 
mission line resonator are found by solving the generalized 
eigenproblem Va„ = w^Ta„ with normalization condition 



ajTa^ = (5 M „ 16011 . In the new coordinates <f> = J2i 4>i e i = 
Cf a ^ tne resonator Lagrangian takes the simple form 



-Y 

2 ^ 



<4® 



(33) 



where v = 0, 1, 2, . . . enumerates the resonator modes. 

In our case, the kinetic matrix T is readily invertible. This 
allows us to further simplify the problem: instead of a general- 
ized eigenproblem, we only need to solve the ordinary eigen- 
value problem 



T^Va, 



w 2 ,a„ 



(34) 



with eigenvector normalization again given by ajTa M = 6^ 
Explicitly, the matrix on the left-hand side of Eq. (f34b reads 

T- X V 



c dz 



c dz 



tc{dzf 



Cl+c dz 



C L +cdz 



(35) 



-1 



V 



2 

c dz 



-i 

c dz 



Cr+cc 



Cr+c dz 



) 



In the continuum limit, where the number of LC elements N 
is sent to infinity and the length of the resonator L = N dz is 
kept constant, the discrete mode vector a„ turns into the con- 
tinuous mode function ip v (z). From the rows i = 2, . . . , (N — 
1) of the matrix equation d34l ), one extracts the second-order 
differential equation 



= -{UvVlc) <p v {z). 



(36) 



The rows i = 1 and i = N yield the homogeneous boundary 
conditions 



dtfiu 
dz 
difv 
dz 



z=0 



= £C l ujI (p v 
= ICrloI Lp v 



z=0 



(37) 



(38) 
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Finally, the orthonormalization condition turns into 

<-L 



z=0 



z=L 



'o 



(39) 



Together, Equations (l36l>-(|39l> form a Sturm-Liouville prob- 
lem 17911 which determines the sinusoidal mode functions 



(40) 



(p v (z) = Acos(uj v v tcz) + B cos(w 1/ V Icz) 



and the corresponding mode frequencies w„. The frequencies 
are obtained as solutions of the transcendental equation 



tanw = 



(xl + xr)u 



1 



XLXR^ 



(41) 



where uj = cjV IcL and \ a = C a /(cL). We emphasize 
that the treatment presented in this section has been exact and 
no assumptions have been made regarding the strength of the 
coupling between the resonator and the left and right circuits. 
In total, the exact Lagrangian (l29l can be written in terms of 
transmission-line eigenmodes as 



a=L.R 



(42) 



B. Model for array of resonators and coupling elements 

For the derivation of the Hamiltonian describing an array 
of resonators coupled by identical superconducting circuits at 
resonator junctions (see Fig. O, we consider the regime of 
weak coupling, as realized in the majority of circuit QED ex- 
periments. Specifically, we will assume that the coupling ca- 
pacitors C c (here, C c = Cl = Cr), connecting transmission- 
line resonators and coupling circuits, are small compared to 
the total capacitance of the resonator, i.e. C c <C cL. In this 
weak-coupling regime, the Hamiltonian takes a particularly 
simple and intuitive form as we shall demonstrate in the fol- 
lowing. 

Quite generally, the Lagrangian of the array can be written 

as 



(43) 



where the terms describe the transmission-line resonators 
("tl"), the ring circuits embedded in the resonator junc- 
tions ("ri"), and the interaction between them ("int"), respec- 
tively. As shown in the previous subsection, the resonator La- 
grangian can be written in terms of eigenmodes v = 0, 1, . . . 
as 



tl. A 



(44) 



We note that for small ratios C c / (cL) the transcendental equa- 
tion fiTI ) can be solved approximately, and the lowest modes 




FIG. 9: (Color online) (a) Array consisting of transmission-line res- 
onators and coupling circuits in the junctions between resonators. 
The coupling circuits, attached to the resonators by capacitors C c , 
are Josephson rings, see panel (b). They consist of a superconducting 
ring interrupted by three identical Josephson junctions with Joseph- 
son energy Ej and junction capacitance C.j. By applying an external 
magnetic field perpendicular to the plane, the loops may additionally 
be threaded by a magnetic flux $. 



are given by cj„ rj vlu - Here, the fundamental frequency 
corresponds to the A/2 resonance and is given by lj /2tt = 

The coupling elements, which will be realized as small su- 
perconducting circuits [Fig.|9|b)] and discussed in more detail 
below, have the generic Lagrangian 



1 



£,„. -«/>, Cf, 1,:^.<I.) : 



(45) 



where C is the circuit's capacitance matrix and V collects all 
inductive contributions of the circuit, including the effect of a 
magnetic flux $ applied to the rings. Finally, the capacitive 
interaction between coupling circuits and resonators is given 
by 



£int,Aj = - m Aj C c(e^ Aj <^-)E^^( ZA i)- 



(46) 



Here, rrtAj plays the role of an adjacency matrix which con- 
tains all information about which resonators are coupled to 
which rings. It is hence defined as 



1 if resonator A couples to ring j, 
otherwise. 



(47) 



Since each ring consists of three superconducting islands, we 
further define a component function fixj £ {1,2,3} which 
selects the individual degree of freedom involved in the cou- 
pling between ring j and resonator A; e MA is the correspond- 
ing three-component unit vector. The coupling capacitors 
(assumed identical across the array) are denoted by C c , and 
z\j = 0, L gives the z variable entering the resonator mode 
function Lp v [as defined in the previous subsection, Eqs. 
©]. 
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To put the circuit and resonator variables on equal foot- 
ing, it is convenient to temporarily rescale the circuit variables 



4>, 



_]/2 • 

C F a so that C,\ v and Fj have identical dimensions. 



C has dimensions of a capacitance, and its magnitude is cho- 
sen such that the nonzero entries in the rescaled capacitance 
matrix K a = C~ 1 Q a are of order unity. 

With these preparations it is possible to obtain an approxi- 
mate expression for the Hamiltonian describing the resonator 
array coupled via Josephson rings. First, the conjugate mo- 
menta are obtained as 



dC C c . T 



dC 



a 



KFj-^mAj^^e^Ou^^Aj). (48) 

*>u 



X.,' 



The coupling terms on the right-hand side of the last two equa- 
tions are small in the weak-coupling limit, C c j \JC Q cL <C 1 
valid whenever cL C c , C D [80] . The inverse of Eqs. (08]), 
required for the Legendre transform, can then be approxi- 
mated by 



Cai/ ~ IXv + 



x c 

m Aj ■-?= ( e u A , K " 1 Qj ) <Pv ( z Xj ) , 



3 



X,v 



C c 

s/C~ 



(K-'e^Kx^M- (49) 



C. Josephson ring couplers 

The coupling elements [see Fig.|9jb)] are located in the res- 
onator junctions and are composed of superconducting loops, 
each interrupted by three identical Josephson junctions. By 
applying an external magnetic field B, each loop may be 
threaded by a magnetic flux $. For reasons to be detailed 
below, we additionally consider the possibility of tuning the 
electric potential of the three superconducting islands by cou- 
pling them capacitively (C g ) to gate voltage sources. The 
Hamiltonian for one such coupling circuit is then given by 



c b -) T C- 1 (Q i -q 3 -) + ^ ! $), (54) 



where the charge vector Qj = [(),. : . (J ,. >■ Qj.3) collects the 
charges on nodes // = 1,2, and 3 of Josephson ring number j. 
Similarly, = C g Vj is composed of the corresponding off- 
set charges. The first term thus represents the ring's charging 
energy and involves the inverse of the capacitance matrix 



C = 



Ce —Cj 
-Cj Ce 

-Cj -Cj 



(55) 



built from the junction capacitances Cj and the sum capaci- 
tances Cs = 2Cj + C c + C g . The inductive energy contribu- 
tions are given by 



In these last equations, we have retained the leading order, and 
corrections are of the order of 0(C^/[C cL]). As a result, the 
weak-coupling Hamiltonian can be written in the form 



H — ^ Ha,x 
x 



E 

3 



Hri.j + ^ ' -Hint, A j 7 



(50) 



with 

#tl,A 

and 



-y 



y^v{a\ v a Xv + -) (51) 



(52) 



(Note that we have reverted back from our temporary rescal 

— 1 /2 - 

ing and that Q j = C Qj has proper dimensions of electrii 
charge.) Finally, the coupling Hamiltonian is given by 



Hmt,\,j = mxjC c (e^ Xj C 1 Qj)yq\„(Pv(z\j), (53) 

The form of the coupling Hamiltonian obtained with Eq. (l53l 
has a simple interpretation: the voltage Vj M = eJ^.C^Qj 
of coupling element j (component /1) is coupled by the ca- 
pacitor C c to the voltage ~^2 V qxv<£v{zxj) at the correspond- 
ing end of resonator A. It is important to note that in the 
Hamiltonian formalism, this intuitive form of the coupling is 
valid only in the weak-coupling limit. As soon as higher- 
order terms 0(C^/[C cL]) are included, the coupling be- 
comes more complicated. 



-Ej Y~] cos 



2n 

— (^ >+1 -^ > -$/3) 

1*0 



(56) 



where the \x indices, enumerating the superconducting islands 
within one ring j, are understood modulo 3, i.e. fi + 1 = 4 and 
\x = 1 are to be identified. For the following discussion, it is 
convenient to drop the ring index "j" and to switch to dimen- 
sionless charge and flux variables defined by = Q p /(2e), 
= 27T0 At /$ o , and tp = 2tt$/$ - 

It is intuitively clear that the total charge N = ni + ni + ri3 
on each ring is a conserved quantity. Formally, this can 
be confirmed by demonstrating that the total charge opera- 
tor and the ring Hamiltonian commute, i.e., using the canoni- 
cal commutators [n^, e ±tLp ^'] — T^W' 6 ^^ one verifies that 
[N, H„] = holds. The eigenstates of the Josephson ring 
Hamiltonian can consequently be written in the form \N,k), 
where k = 0, 1, . . . enumerates the eigenstates in the subspace 
of total charge N. 

We assume that a residual coupling of the circuit to its en- 
vironment allows it to relax into its ground state \ipo) = 
I No , ) . Noting that the interaction Hamiltonian H mt also 
commutes with N, we will assume that, for the duration of 
an experiment, the circuit remains in this ground state. The 
virtual intermediate states involved in the transfer of photons 
correspondingly belong to the same total charge subspace and 
hence can be written as | Nq, k ). 

Since, in the general case, the ring Hamiltonian is not 
amenable to an analytical solution, we obtain its spectrum and 
the charge matrix elements (required in the subsequent sub- 
section) by numerically exact diagonalization. Our strategy 
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FIG. 10: (Color online) Dependence of the ground state charge num- 
ber Ao on external magnetic flux $ and offset charges, here for the 
uniform case n s i = n 9 2 = n g 3 = n g . As expected, No takes 
on only integer values corresponding to the total number of extra 
Cooper pairs located on the Josephson ring. The integer-step bound- 
aries between regions of different Ao in general acquire a small fi- 
nite width due to the residual coupling to the environment that allows 
charge relaxation. Parameters chosen for this plot: Ej/h = 10 GHz, 
C.j = 0.7 fF, and C c = 5fF, yielding Ej/E s ~ 2. 



is as follows: in the first step, we employ diagonalization in 
the charge basis to obtain the ground state | i[>o ) and use it to 
extract the total charge, 

iV = (Vo|iV|Vo). (57) 

Numerical results for this ground state charge in a Josephson 
ring with realistic parameters are presented in Fig. |T0j As 
can be inferred from the figure, Ao is generally an integer- 
valued function of both offset charges and external magnetic 
flux. In the regime of strong charging effects, the dependence 
on flux weakens, and explicit expressions can be obtained for 
the boundaries between No regions in offset-charge space [see 
Appendixlcl. 

In the second step, we may then restrict ourselves to one 
particular subspace of total charge No. To do so, we perform 
a canonical transformation 

If! = (f[ + <p' 3 , lfi 2 = tp' 3 - f2, = </?3 ( 58 ) 

7ii = n[, n 2 = —n' 2 , n 3 = -n[ + n' 2 + n' 3 , (59) 

after which the variable cp' 3 is cyclic and the corresponding 
canonical momentum ng = n\+ri2+n 3 = N is the conserved 
total charge. With this, the restriction of the Hamiltonian to 
the Ao subspace can be brought into the form 



H^ No) =4Ex ( n[ - ±[n fl i - n g3 + 



lg2 



lg3 



N 



(60) 



Ej cos ( y?i - ^ j - Ej cos ( tp' 2 - | j - Ej cos f ip' x + tp' 2 + | 



Here, the charging energy Ey, has been defined such that 
4i?s = (2e) 2 (7i — 72), and 71.2 are reciprocal capacitances 
obtained in the inversion of the capacitance matrix C, see Ap- 
pendix [B] The Hamiltonian H^ ^ has one degree of freedom 
less than the original ring Hamiltonian H r i, and is thus more 
convenient for the numerical calculation of eigenenergies and 
charge matrix elements. 

In preparation for the next subsection where the Josephson 
rings will be integrated out (relying on the dispersive limit), 
we finally rewrite the interaction Hamiltonian in the subspace 
Ao- For the example of a single Josephson ring coupled to 
three resonators, the component function [i in the coupling 
Hamiltonian (l53l l takes the simple form = j Con- 
sidering only one of the low-lying modes of the resonators, we 
will drop the mode index "v" from here on, and write u) r for 
the (angular) resonance frequency. For the coupling Hamilto- 
nian we then obtain 

flint = C c V nm {a + at) T C- 1 Q, (61) 

where the vector a collects the annihilators for the three res- 
onators A = 1,2,3, which are obtained by rewriting q\ = 



y / oj r /2(a\ + a\). Vans = yJuj r /2Lp{0) w u r / cL is the 
root-mean-sguare voltage in the resonators at the relevant res- 
onator end (8lj|. Once the Hamiltonian doTT ) is restricted to the 
subspace of total charge Ao, one can show that it assumes the 
form 

H^ o) =2epV ms n' 1 {a 1 - a 3 + H.c.) 

+2epV rms n' 2 (a 3 - a 2 + H.c.) (62) 

with capacitance ratio /3 = C c (7i — 72). Note that here we 
have discarded terms of the form a(a\ + at) with a repre- 
senting a c-number. Such terms merely displace the resonator 
mode, and can ultimately be absorbed into a redefinition of 
the offset charges. 



IV. EFFECTIVE PHOTON LATTICE HAMILTONIAN 

We now turn to the crucial step of integrating out the 
Josephson ring elements and specifying the conditions under 
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which the resulting photon lattice Hamiltonian breaks time re- 
versal symmetry. The adiabatic elimination of the degrees of 
freedom of the coupling circuits is based on being in the dis- 
persive regime of large energy mismatch between photonic 
excitations of the resonators, and excitations of the coupling 
circuits. Specifically, the dispersive regime is defined by the 
inequality g <C A, where A represents the detuning between 
photonic and circuit excitations and g is the effective strength 
of their mutual coupling. For a general and systematic exposi- 
tion of the adiabatic elimination technique we refer the reader 
toRef. 13^1 . 

Working within the rotating-wave approximation (RWA), 
the total number of (dressed) photons is conserved. For a 
given total photon number, we define Pq as the projector 
(Pq = 1) onto the subspace with that photon number and 
with all Josephson rings occupying their ground states. The 
effective photon lattice Hamiltonian H v ^ can be obtained by a 
canonical transformation, 

H ph = P e lS He- lS P (63) 
= ]T H ±x + ^P [iS, H mt ]P + 0{Hl t ), 

A 

where the generator S of the transformation is chosen such 
that the linear coupling between rings and resonators is elimi- 
nated. To leading order in the interaction, it is given by 

^E ^ ^la'Xaifi-Hx. (64) 

a, a' 

where a, a' are indices for the eigenstates of H t \ + in the 
P subspace, and Pi = 1 — P projects onto the comple- 
mentary subspace. The main task hence consists of evaluating 
the contribution hPo[iS, iJ; nt ]Po to the effective Hamiltonian. 
Following the arguments about charge relaxation in the pre- 
vious subsection, we carry out this evaluation in the subspace 
with charge No, which contains the ground state of the cou- 
pling elements. 

To illustrate our procedure, we consider the simple case of 
three resonators attached to a single coupling element. [The 
generalization to a full array can be achieved by starting from 
Eq. d53l and projecting it onto the No charge subspace of all 
rings.] In RWA, the interaction Hamiltonian d62t reads 



fe>0 



ni,jk|JVo,fc)(JV 0j O|(ai-a3) 



+ n 2 , k \N Q ,k)(N ,0\(a 3 -a 2 ) 



He, (65) 



where n^k — ( No, k \ n' \ No, ) denotes the relevant charge 
matrix element. It is crucial to note that the origin of photon 
hopping with complex-valued hopping elements is directly 
based on the fact that these charge matrix elements may be 
non-real, as we will see momentarily. A tedious but elemen- 
tary calculation shows that the effective photon Hamiltonian 
is given by 



H ph = y^(a; r + e\)a\a\ + 



A=l 



3 

E 

A=l 



i A a A aA+i +H.c. 



(66) 



where the index A in the second term is to be understood as 
A mod 3, and where the energy shifts and photon hopping ma- 
trix elements are found to be 



ei = 2(/3^ ms ) 2 £^l 
^ LJ r - E k 



k>0 



e 2 = 2{peV rms ) 2 J2 



2V l n 2,fc| 5 



k>0 



(67) 
(68) 



e 3 = 2WeV lms ?Y, lni ' k ~T l \ 



k>0 



t 1 =2(/3eV lms ) 2 Y, 



fc>0 



ta = 2(PeV uns ) 2 



fc>0 



LU r — Ek 

-(ni.k)*n 2 ,k 
io r — Ek 

w r — Ek 



(70) 



, (71) 



* Of«„T/ \2 ( n l.k)* n 2,k ~ 1^1, A:| 2 



k>0 



ui r — Ek 



Ek denotes the energy of the k-th circuit excitation (measured 
relative to the ground state energy Eq). Eqs. d70l>— d72b for 
the hopping matrix elements confirm our previous statement 
that the emergence of complex phase factors in the hopping is 
directly linked to the possibility of non-real charge matrix ele- 
ments. Before investigating the conditions under which these 
charge matrix elements are non-real and result in breaking of 
time-reversal symmetry, it is useful to note that, in general the 
above equations will also lead to breaking of the three-fold ro- 
tation symmetry due to the energy shifts e\. The origin of this 
is, of course, the possible presence of different offset charges 
on each of the three superconducting islands. 

For the present discussion, we restrict our discussion to 
the case where no such breaking of the three-fold symme- 
try occurs, and we will hence choose identical offset charges 
rigi = n g2 = n g 3 = n g . In the ideal case, individual super- 
conducting islands would not need to be connected to separate 
gate voltage sources; instead, a global electric field perpendic- 
ular to the chip plane (e.g., by a back gate) could be applied to 
achieve a uniform and tunable offset charge. (This, of course, 
neglects the presence of random offset charges and 1 // charge 
noise which we address in Section [IV CI ) With the threefold 
symmetry intact, one concludes that 



(l = £2 = e.3 



(73) 



must be satisfied. In other words, application of a global elec- 
tric field does not lead to energy detuning between resonators. 

We need to be cautious though not to throw out the baby 
with the bath water. Clearly, fixing all offset charges to be 
identical is a strong restriction of parameter space and it is by 
no means obvious that this leaves any freedom for complex- 
valued matrix elements and hence time-reversal symmetry 
breaking on the level of the effective photon Hamiltonian. Let 
us thus verify that Eq. ( f73l when combined with Eqs. d67ii- 
(|72| i is in general compatible with complex-valued hopping el- 
ements t\. Given that e\ must take the form of Eqs. <f67l>— (l69b. 
a sufficient condition for satisfying €\ = e 2 = £3 is obtained 
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by requiring that, for each excitation level k, the charge ma- 
trix elements n^k have equal modulus, \n\ k\ = l n 2,fc|> and 
obey |ni,fc| 2 = |ni,fc — n2,k\ ■ Evaluating these conditions, 
we find that the charge matrix elements obey 

n^ k = r k e if »' h (74) 

with modulus ru > independent of the charge index // = 
1, 2, and phases 

/l,*-/2,* = (±)*|+27TZ fc . (75) 

The latter equation must hold for all levels k = 1,2,..., but 
both the sign and the integer € Z may differ among lev- 
els. The freedom in the phase sign turns out to be crucial 
for breaking time-reversal symmetry. Without the sign free- 
dom or when truncating the system to a two-level system, 
the (gauge-invariant) phase sum over the three-resonator loop 
would always be an integer multiple of tt. Hence, as discussed 
in Section []TJ time-reversal symmetry would be intact on the 
level of the effective photon Hamiltonian. However, due to 
sign flips for higher levels k and together with the different 
prefactors in the terms of the sum [Eqs. (|7Q]i-(|72]il, arbitrary 
gauge-invariant phases 

3 3 

£ = arg Y[ t\ (76) 

A=l A=l 

can in principle be generated and time-reversal symmetry thus 
be broken. 



A. Numerical results for intermediate E j /Ey. 

Equations (l70ii-(l72ii allow for a direct evaluation of the 
essential parameters of the effective photon Hamiltonian. 
The most important quantity for determining whether time- 
reversal symmetry breaking succeeds is the gauge-invariant 
phase sum^ c <p, Eq. (|76"1 >. Whenever this sum corresponds to 
an integer multiple of it, time-reversal symmetry is intact; for 
all other values it is broken. In these terms, our prime concern 
is to demonstrate that 

§T <p i ttZ (77) 
c 

can be achieved for realistic device parameters and reason- 
able magnitude of the photon hopping element (clearly, for 
hopping matrix elements with \t\ = the complex phase be- 
comes arbitrary and completely meaningless). 

Results from numerical diagonalization for a selected set of 
parameters, chosen with current fabrication capabilities and 
general parameter requirements in mind, are presented in Fig. 
[TT1 We find that breaking time-reversal symmetry is feasible 
under realistic conditions, and that the external dc electric and 
magnetic fields can be utilized to switch time-reversal invari- 
ance on and off (with the electric field) and to smoothly tune 
the value of the gauge-invariant phase sum (with the magnetic 
field). 



Several comments are in order to provide an intuitive under- 
standing of the numerical results shown in Fig. [TT] We note 
that the excitation energies of the Josephson ring and the re- 
sulting photon hopping amplitudes and phases exhibit a step- 
like dependence on the global offset charge. This is easily un- 
derstood from the Josephson ring Hamiltonian, Eq. d60l i: The 
values of the offset charges fix the total charge Nq. Further, 
in the case of identical offset charges n g \ = n g 2 = n g s, this 
is the only way the offset charges enter the Hamiltonian. By 
consequence, the fact that Nq is an integer-valued function 
of n gfi explains the step-wise dependence on offset charges. 
Only at special points where an increase in the common offset 
charge causes a level crossing of the two lowest states in sub- 
spaces with different total charge, the parameter No changes 
discontinuously from one integer to another and thus leads to 
the observed steps. 

The fact that time-reversal symmetry is broken for No = 
1 , 2 (and, by means of charge periodicity, for all No mod 3 = 
±1) and that the gauge-invariant phase sums are of opposite 
sign for these two cases can easily be motivated by consider- 
ing the case of large charging energy. For Nq = 1 there are 
three nearly degenerate states with one additional Cooper pair 
(the "particle") located on one of the three islands. When Ej 
is finite, the extra Cooper pair can start to move, becomes sus- 
ceptible to the vector potential and produces an effective phase 
in the photon hopping. Conversely, for No = —1 (equivalent 
to No = 2) there are three nearly degenerate states with a 
Cooper pair missing (i.e., a "hole") on one of the three is- 
lands. This results in the opposite signs of the gauge-invariant 
phase sums since hopping of particles involves the phase <p, 
whereas hopping of holes is associated with phase —(p. The 
case Nq mod 3 = corresponds to the particle-hole sym- 
metric case, where the photons acquire zero synthetic gauge 
charge and time-reversal symmetry holds. 

As we will prove below, the regime of very large Ej / 
ratios (where Josephson tunneling completely overwhelms 
charging effects) is inadequate for breaking time-reversal 
symmetry. As a result, charge noise must be expected to im- 
pose limitations on the proposed device, which we briefly ad- 
dress in Section lTV CI Future work must establish the optimal 
working point where ip comfortably reaches the crucial 
value of 3 x ir/6 = ir/2 while keeping sensitivity to offset- 
charge fluctuations at a minimum. 



B. Conditions for time-reversal symmetry breaking 

First, let us establish that in the regime where Josephson 
tunneling dominates over charging effects, i.e., Ej/E^ ^> 1, 
the Josephson ring fails to break time-reversal symmetry. To 
see this, consider the ring Hamiltonian d60b in phase basis 
where n! = id/ dip' (we will drop primes in the following). 
For Ej 3> Ex, the Hamiltonian describes the situation of a 
fictitious particle with large mass in a two-dimensional poten- 
tial. (Strictly speaking, the space described by the coordinates 
tpip is a torus, since the periodic boundary conditions require 
that ipp and ip^ + 2ir be identified as the same coordinate.) 
Due to the large mass, the low-energy part of the spectrum can 
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FIG. 11: (Color online) Numerical results for a junction of three resonators attached to a central Josephson ring coupler. The device is tunable 
by varying the magnetic flux see color/gray scale, and by changing the global offset charge n g as set by a constant electric field, see 
x axes. Panel (a) shows the lowest transition frequency o>oi/27r of the Josephson ring device in comparison with the resonator frequency 
u} r /2n — 7 GHz. As one can check, the dispersive limit is maintained for the selected values of magnetic flux $. Panel (b) displays the 
resulting magnitude of photon hopping strengths \t\. The non-monotonic behavior is explained by the crossing of the ujoi transition and the 
resonator frequency around <&/$o ~ 0.3. Panel (c) presents the corresponding results for the gauge-invariant phase sum^ tp and proves the 
breaking of time-reversal symmetry. As expected from general considerations, time-reversal invariance remains intact at zero offset charge, 
and at zero magnetic flux. (Parameters as in Fig.[lO] in addition: C c = 5 fF, oj,./2ty = 7 GHz and \/T[c = 50 fi.) 



be described by a local approximation of the two-dimensional 
potential at its minimum [82], 

V(<p) ~ i(v> - <p min ) T M{f - (prnrn)- (78) 

Here, M is positive definite, and we have used the vector no- 
tation ip = ((fi,(fi2)- (Note that both the curvature matrix 
M and the position of the minimum tp m i n still depend on the 
magnetic flux, which we suppress in our notation.) Once 
the approximation (l78l l is employed, the periodic boundary 
conditions are changed into the regular boundary condition 
f R2 dtpidip2 1^0(^1 , V?2)| 2 = 1- This opens the way for a 
gauge transformation 

ip(<Pi,<P2) = exp(iai^i + ia 2 <f2)'4>(<Pi, (79) 

which leaves the new boundary condition unchanged. Choos- 
ing 

a m = {-l) m {N + 3n gm - ££ =1 %M )/3, (80) 

this transformation can be used to eliminate all offset-charge 
related first derivatives from the Schrodinger equation for ■ip. 
In other words, in this gauge the fictitious particle does not 
"see" a vector potential and its ip wavefunction can be cho- 
sen entirely real-valued. This in turn reveals that all charge 
matrix elements can be chosen purely imaginary, and conse- 
quently all hopping elements for photons purely real-valued, 
tfj. £ R JHI. While time-reversal symmetry is thus not broken 
in this regime, we emphasize that Josephson rings in the large 
Ej/Ey, regime are still very useful: they make the photon 
hopping strength tunable with an external magnetic field 
and remain insensitive to the effects of random offset charges 
and 1// charge noise just like the transmon qubit I6lll62ll . 

Closely related to the no-go statement for time-reversal 
symmetry breaking with large Ej / E-g ratios, one can specify 



two general conditions required for breaking of time-reversal 
symmetry. First, we note that breaking particle-hole symme- 
try, or equivalently, the presence of nonzero offset charges, is 
required. The argument for this directly follows from our pre- 
vious discussion: without offset charges, all eigenfunctions 
of the Josephson ring Hamiltonian in phase basis can be cho- 
sen real-valued outright [i.e., without the substep of approxi- 
mating the potential in Eq. d78lll. The repetition of our argu- 
ments following Eq. (l80l then again leads to the conclusion 
of no time-reversal symmetry breaking. For the case of iden- 
tical offset charges, we can narrow down the necessary con- 
dition further: since the Hamiltonian (|60l remains invariant 
(up to an irrelevant overall constant) under the transformation 
jVo —> Nq ± 3, we find that Nq mod 3 = ±1 is required to 
break time-reversal symmetry. 

Second, we note that the presence of Josephson junctions 
is crucial in our scheme. Without them, the inductive energy 
would generically take the form of Eq. ( l78l . and all subse- 
quent arguments leading to the conclusion of no time-reversal 
symmetry breaking hold. 



C. Consequences of random offset charges and 1 // charge 
noise 

It is known from experiments with superconducting charge 
qubits I63M66I1 that the coupling of a superconducting circuit 
to its environment generally results in random offset charges 
on superconducting islands, and that these offset charges typ- 
ically fluctuate as a function of time with a characteristic 1// 
noise spectrum ll66l - l68ll . This behavior will likely affect the 
performance of the Josephson coupler circuits proposed here, 
and we comment on consequences and potential solutions to 
this issue. 

For superconducting charge qubits, the negative effects of 
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FIG. 12: (Color online) (a) Josephson ring with attached voltage 
bias lines for cancelling random offset charges, (b) Effective pho- 
ton hopping strengths and gauge-invariant phase sums for random 
offset charges, with n 9M £ [0, 1] with uniform probability distribu- 
tion. Data points are placed such that their x positions correspond 
to the gauge-invariant phase sums ip (modulo 2ir), and their y 
positions display the arithmetic mean of the three photon hopping 
strengths |t M |. For each data point, an "error" bar shows the spread 
from the minimum |t M | to the maximum. (Device parameters used 
are the same as in Fig.ll II) 



charge noise can be cured by working with transmon qubits in 
the regime where Josephson tunneling dominates over charg- 
ing effects I6lll62ll . This venue, however, is not available for 
the Josephson ring circuit when aiming at time-reversal sym- 
metry breaking, as follows from our discussion in the previous 
section. While devices with large Ej / Ey, will be insensitive 
to charge noise and very useful for making photon hopping 
strengths tunable, the gauge-invariant phase sum around the 
loop will be exponentially suppressed. 

For devices with one or maximally a few Josephson coupler 
circuits, it is conceivable to work with intermediate Ej /Ey 
ratios and to couple the individual superconducting islands 
capacitively to voltage bias lines, see Fig. [12 a). This way, 
random offset charges can be cancelled and the device sta- 
bilized. For larger arrays, attaching individual bias lines be- 
comes cumbersome. Random offset charges then lead to dis- 
order in the photon hopping elements as well as in the gauge- 
invariant phase sums, see Fig. [T2l b). While presence of such 
disorder poses interesting questions itself (compare, the re- 
cent interest in potential disorder in ultracold atom systems, 
see e.g. ll69tl . and localization in random magnetic fields, e.g. 
[TtoL l7llo . future studies will also aim at identifying alterna- 
tive superconducting circuits for charge-noise insensitive and 
time-reversal symmetry breaking coupling elements. 



V. CONCLUSIONS AND OUTLOOK 

In summary, we have shown that superconducting circuits 
based on Josephson junctions can be used to break time- 
reversal symmetry in arrays of on-chip microwave resonators. 
In the first part of our paper, we have explored how to use pas- 
sive coupling elements to generate gauge-invariant phases in 



the lattice hopping elements, and how these phases are related 
to time-reversal symmetry breaking. Much of this discussion 
is general and can readily be transferred to lattices other than 
photon lattices. Our subsequent discussion has highlighted 
consequences and applications of breaking time-reversal sym- 
metry in non-interacting lattices of photons, including the re- 
alization of an on-chip circulator and the achievement of a 
highly tunable band structure for the concrete case of a pho- 
tonic Kagome lattice. We note that the existence of localized 
photon states on hexagons in the Kagome lattice may be of 
interest for photon storage in the future. These localized pho- 
ton states do not necessitate the presence of a large lattice, but 
can rather be accessed in a single Kagome star consisting of 
only twelve resonators - a setting that is well within reach of 
current experimental capabilities. 

The second part of our paper has addressed a concrete pro- 
posal for the realization of such passive coupling elements 
in the circuit QED architecture. Our presentation aimed to 
be pedagogical and to collect the necessary circuit quantiza- 
tion tools to handle an array of transmission line resonators 
coupled to small superconducting circuits playing the role of 
coupling elements. We have stated the general conditions for 
breaking time-reversal symmetry with a passive coupling ele- 
ment, including the necessity of non-linear elements (Joseph- 
son junctions), the presence of a magnetic field, and break- 
ing of particle-hole symmetry. We have shown that an ex- 
tremely simple circuit, a superconducting ring interrupted by 
three Josephson junctions, can be used to satisfy all the nec- 
essary requirements. For realistic device parameters, we have 
calculated the resulting photon hopping strengths and gauge- 
invariant phases as a function of external magnetic flux and 
global offset charge. Finally, we have identified random off- 
set charges and charge noise as likely challenges when target- 
ing a lattice without disorder in hopping strengths and phases. 
Future works will explore alternative circuits for tackling this 
issue, and will address the interesting question of strongly- 
correlated photon states with broken time-reversal symmetry, 
which are expected for large effective photon-photon interac- 
tion such as in the Jaynes-Cummings lattice. 
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Appendix A: Time-reversal symmetry 

Generally, the dynamics of a system is said to be time- 
reversal symmetric if for a given solution to the equations 
of motion, the corresponding motion-reversed evolution is a 
valid solution as well. In the following, we briefly compile 
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the most important facts about time-reversal in quantum me- 
chanics. 

In quantum mechanics, symmetries manifest as maps S of 
Hilbert space, which leave all observable probabilities invari- 
ant i.e. | ( | Sip ) | 2 = |(</>| V)| 2 for all states \(f>) , \ tp) 
[F72L F73h . This is fulfilled if and only if S is either a unitary 
operator, or an operator which is anti-linear and anti-unitary 
JzUll- While the for mer choice applies to discrete and con- 
tinuous symmetries including rotations and parity, the latter 
option must be selected for time reversal, in order to avoid en- 
ergy spectra not bounded from below [see, e.g., Ref.|72|for the 
proof of this statement]. The time-reversal operation must 
thus be anti-linear and anti-unitary, i.e. 



@(a\(f>) +(3\ V>) =a*6 

(e0|evo = <vi0>- 



0> + /3*e| 



(Al) 
(A2) 



Once time reversal 9 has been properly defined for a specific 
system with Hamiltonian H , symmetry of that system under 
time reversal is signalled by the fact that QHQ^ 1 = H holds. 
(For simplicity, we are excluding the case of degenerate eigen- 
states of H, for which may additionally induce a rotation 
within the degenerate subspace.) 

To define explicitly, we assume that the system provides 
us with an observable (with non-degenerate spectrum), say x, 
which is expected to be time-reversal invariant for physical 
reasons. For example, this operator may be the position op- 
erator for the location of a particle in real space; for a circuit 
network, it may be the operator for charge on a certain net- 
work node, which also must remain invariant under time re- 
versal. Under these assumptions, time-reversal is expected to 
leave the eigenstates of x invariant, possibly up to a phase, 







(A3) 



from which 0x0 _1 = x immediately follows. Time-reversal 
symmetry thus holds if and only if there exists a phase $(x) 
such that 0iJ0 _1 = H is satisfied. We will see momentarily 
that the phase 1? is intimately related to phases arising from 
gauge transformations. 

Eq. dA3t has several important consequences, which we 
briefly gather in the following, (i) Once $(x) is fixed, the 
action of on the entire Hilbert space is uniquely defined by 
Eq. jA3\ . To see this, decompose any state | ip ) in the position 
basis and invoke anti-linearity to obtain 



6|V) 



d d xQ 



V>(x)|x) 



d d xe 4,5(x V*(x)|x) 



(A4) 



(ii) The anti-unitarity condition, Eq. (IA21) . is automatically 
satisfied by this definition of 0. (iii) The canonical momen- 
tum p transforms under time-reversal as 



0P0- 1 



W(x), 



(A5) 



which can be derived using Eq. dA3t and the canonical com- 
mutator [x, p] = i. 



To demonstrate how the phase $ is determined by our gauge 
choice, consider the example of a particle with mass m in an 
external potential with Hamiltonian H = p 2 /2m + V(x). 
Choosing i?(x) = 0, one can verify that 0p0 _1 = p, 
and hence 0i70 _1 = H. As expected, the problem is time- 
reversal symmetric. The same system can, of course, be de- 
scribed in a different basis, related to the original position ba- 
sis by a local gauge transformation, | x ) n- e* x ( x ' | x ) . In the 
transformed basis, the Hamiltonian takes the modified form 



H 



1 

2m 



P + Vx(x) 



(A6) 



Performing a gauge transformation cannot affect time-reversal 
invariance, and so 0i70 _1 = H should hold for an appropri- 
ate choice of fl. Indeed, using Eq. 1A51 . we can construct -d 
by requiring 



H = QHQ- 1 = 



1 

2m 



p + W(x) + Vx(x) 



V(x) 



(A7) 

which yields V^(x) + 2V%(x) =0. As a result, the phase 
of the time-reversal operator is fixed by the gauge, ??(x) = 
— 2%(x) up to an irrelevant constant. If we interpret A = Vx 
as a vector potential (here with zero curl), we can write 



tf(x) = -2 ds ■ A 



(A8) 



As an immediate corollary we note that the presence of a 
magnetic field would manifest in a vector potential A with 
nonzero curl. In that case, the resulting equation Vi?(x) + 
2A = has no solutions, and hence time-reversal symmetry 
is broken. 

In summary, one can thus show that the following equiv- 
alences hold for the case of position and momentum oper- 
ator having continuous spectra: Time-reversal symmetry is 
intact. There exists a phase choice for i9(x) such that 
0i70 _1 = H holds. There exists a local gauge trans- 
formation that makes the Hamiltonian real-valued. The 
vector potential satisfies § c ds ■ A = for any closed contour 
C. (Note that non-singularity of the phase functions is implied 
everywhere.) 

Finally, let us switch to the case of a discrete position oper- 
ator, such as for a lattice Hamiltonian 



H= |t| JV^it 



kj 



(A9) 

describing a system of particles which can hop between lattice 
sites, say from j to k, and doing so pick up a phase factor 
tfjk- As the analog of the continuous position basis, we use 
the particle number states | nx,ri2, . ■ .}, and hence define the 
time-reversal operation via 



e\ ni ,n 2 ,...) =e M ^^-^\ nil n 2l .. 



(A10) 



For our purposes it is sufficient to consider linear functions of 
the form ri2, . . .) = J2j$j n j- Invariance under time 
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continuous 



discrete 



|x) 

0(x) 

A(x) 

Vtf + 2A = 
/ cfe • A = 



| 711,712,...) 

ti k - ^ + 2<#y <E 2tt2 



TABLE II: Correspondences for time-reversal symmetry in continu- 
ous and discrete systems. The statements in the last two rows only 
hold if the system is time-reversal invariant. 



reversal is then equivalent (by definition) to the existence of 
phases §j such that QHQ^ 1 = H holds. 

From Eq. ( 1A10I ) with linear «?, one obtains the transforma- 
tion law for annihilation operators, which reads 



to-i 



(All) 

Applying the time-reversal operation to the Hamiltonian jA9i , 
we thus find that invariance under time-reversal implies the 
existence of a set of phases {flj} such that 



■dj + 2<p kj G 2?rZ 



(A12) 



holds for all lattice indices j, k. (Note: once such phases 
i9j have been found, the gauge transformation with phases 
{•dj/2} makes the number-basis Hamiltonian real-valued.) 
The last condition ( IA12b can finally be shown to be equiva- 
lent to the requirement that 



ip jk e ttZ 



(A13) 



cm 



for all closed loops C. The correspondences between the con- 
tinuous and the discrete case are summarized in Table HH 



Appendix B: Inverse of the capacitance matrix 

For completeness, we provide explicit expressions for the 
inverse of the capacitance matrix C: 
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The reciprocal capacitances 712 > are defined as 

Ce — Cj 



7i 

72 



(Cx-2Cj)(C s + Cj)-' 
Cj 

(Cv-2Cj)(Cv + Cj)- 



(Bl) 

(B2) 
(B3) 



Appendix C: Total charge number of the Josephson ring ground 

state 



In Section lTlICl we noted that the eigenstates of the Joseph- 
son ring Hamiltonian (154-b can naturally be chosen as simulta- 
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FIG. 13: (Color online) Regions of fixed total charge N = 0, ±1 in 
the charge regime, as a function of the three offset charges n s i, n g 2 
and n ff 3 . The shape of the region boundaries depends on the charging 
energy ratio Ec/E' , chosen as (a) 5/4, (b) 5/2, and (c) 20. Note 
that the coordinate axes are oriented differently in panel (c) to reveal 
the flatness of the boundaries for large Ec/E' c . 



neous eigenstates of the total ring charge N = n\ + n 2 + 713, 
here measured in units of the Cooper pair charge (2e). For 
the subsequent discussion in that section, it was important 
to extract the total charge number of the ground state No = 
(ip \N\tpo) for given offset charges n ff = (n gl , n g2 , n g3 ) 
and model parameters. While numerical diagonalization of 
the Hamiltonian d54l i allows the direct calculation of No, it is 
useful to first understand the general structure of Nq. 

Our starting point is the ring Hamiltonian, written in terms 
of dimensionless charge numbers n = (m, n 2 , H.3) and phase 
difference ipj , 



Hi =4£ c (n-n s ) T M(n-n g ) 



(CI) 



- Ej ^2 cos(ipj - (fj-i - <p/3) = H c + Hj, 

3 = 1 

where AEq — \{2e) 2 l\ is the charging energy associated 
with the reciprocal capacitance 71 [see Eq. ( IB2H . M is a di- 
mensionless matrix obtained from the inverse capacitance ma- 
trix C _1 by rescaling and is defined as (M)jj = (1 — 7)<5jj+7, 
with 7 = 72/71 = E' c /E c - 

Since A" has a discrete spectrum (comprised of all inte- 
gers Z), it is clear that the offset-charge space spanned by 
(n g i , n g 2 , n g z ) is divided into regions of constant ground state 
charge number No- At the boundaries of these regions, No 
must jump discontinuously. To understand the boundaries be- 
tween such regions, we make an important observation which 
is not limited to the charging regime, but holds for arbi- 
trary Ej /Ec ratio, and is also independent of all remaining 
model parameters: Any shift of the offset charges by integer 
amounts, 

rig -> rvg + (zi,z 2 ,z 3 ) (z> e Z) (C2) 



leaves the spectrum of H invariant and shifts Ao according to 

(C3) 



N -> No + z * 



Further, at zero offset charge n g = 0, particle-hole symmetry 
is intact and dictates Ao = 0. From Eq. iC3\ one thus imme- 
diately knows that the ground state charge number obeys 



No(zi,Z 2 ,Z 3 ) = Zx + Z2 + Z3 {z% G Z). 



(C4) 
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Equation iC2i , in fact, allows one to restrict the entire dis- 
cussion to the domain n g j E [—1/2,1/2). Symmetry also 
dictates that, assuming the simplest case of a direct transi- 
tion from 7Vo = at n g = to A*o = ±1 at n g = ±e,-, 
the transition must occur at the midpoints. In other words, 
the points ±(1/2,0,0), ±(0,1/2,0) and ±(0,0,1/2) must 
lie on the boundaries separating iVo = from Ao = ±1. 
Analogous arguments apply for the transition to Nq = ±1 at 
rig = ±ei+e2 — e3 etc. along six out of the eight space diago- 
nals, which puts the corresponding six corners of the unit cube 
on the boundaries. This sets the overall structure of A . The 
detailed form of the full boundary, however, depends on de- 
tails such as the Ej / Ec ratio. In the charge limit (Ej < Ec), 
Nq can be constructed analytically and it is instructive to do 
so and to discuss how Ao is modified for increased Josephson 
tunneling. 

In the charging regime, it is primarily the charging contri- 
bution He which determines the boundaries between Nq re- 
gions. To leading order, we hence neglect Josephson tunnel- 
ing (Hj) completely, and the problem becomes similar to the 
question of charge stability in a triple quantum dot |[75ll . The 
eigenstates of He are charge eigenstates | n ) with n € Z 3 
and corresponding eigenenergies E n (n g ). The boundary be- 



tween the No = region centered at n g = and the adjacent 
Ao = ±1 regions reached via the planar diagonals are ob- 
tained by requiring that the respective energies match, 

E (n g ) = E± ej (n g ). (C5) 

This yields six equations of the form 

= lip 2(1 -i)n gj ±2 7 J2 n ak (.7 = 1,2,3), (C6) 

k 

which define planes in the offset charge space. Consistent 
truncation of the planes to the region where Nq = —> ±1 
can occur, yields the full charge boundaries, see Fig. [13] Note 
that in the charge regime, the Ao boundaries do not depend on 
the magnetic flux. 

The presence of Josephson tunneling will generally mod- 
ify the shape of these boundaries, but leave the properties 
derived from general symmetry arguments intact. We expect 
Hj to introduce flux-dependence and to smoothen the sharp- 
edge boundaries [see, e.g., Fig. [T3T a) and (b)], as it hybridizes 
the states | ej ) for j = 1,2,3 and thus turns crossings into 
avoided crossings. 
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